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We calculate the spin-dependent zero-bias conductance Gg-cr' in armchair graphene nanoribbons with hydro¬ 
gen adsorbates employing a DFT-based ab initio transport formalism including spin-orbit interaction. We find 
that the spin-fiip conductance Gg-a can reach the same order of magnitude as the spin-conserving one, Gg-o-^ due 
to exchange-mediated spin scattering. In contrast, the genuine spin-orbit interaction appears to play a secondary 
role, only. 
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I. INTRODUCTION 

In recent years, graphene [1] has been considered as ideal 
spintronics material [2]: Due to the weak spin-orbit interac¬ 
tion (SOI) [3], spin lifetimes of Dirac electrons are expected 
to be long. However, the Hanle precession measurements typ¬ 
ically yield spin diffusion times several orders of magnitude 
below the theoretical predictions [4] . Recently, a quantum in¬ 
terference measurement [5] proposed that intrinsic local mag¬ 
netic moments at defects [6] are the primary cause of spin re¬ 
laxation in graphene, masking any potential effects of the gen¬ 
uine SOI. The efficiency of such a mechanism was confirmed 
by consecutive theoretical work [7]. While another mecha¬ 
nism for spin-flips originating from bias induced orbital mag¬ 
netism has also been identified, recently, its quantitative effect 
still remains to be explored [8]. 

Motivated by the efficiency of exchange induced spin-flips 
in graphene, we study spin-dependent transport in graphene 
nanoribbons (GNRs), i.e. strips of graphene with ultra-thin 
width, from first principles. Our interest in GNRs is closely 
related to their electronic properties: GNRs inherit a weak in¬ 
trinsic SOI [9] and high electron mobility from graphene [10]. 
Moreover, GNRs exhibit gaps that can be tuned with the rib¬ 
bon width [11] and local spins can be generated at zigzag 
edges [12] or defects [6]. These properties make GNRs 
promising materials for applications in spintronics, e.g. for 
quantum computing [13]. 

Spin transport in GNRs will be addressed in this paper us¬ 
ing the standard formalism of molecular electronics [14]: A 
device scattering region is located between two semi-infinite 
leads with applied bias V\ the total conductance G = dl/dV of 
the device is split into four spin-dependent conductance coef¬ 
ficients Gcrcr' (with G = Yjo-o-' Gcrcr')- They account for an elec¬ 
tronic current being injected with spin cr, and, after passing 
the device region, measured with spin direction cr'. 
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In ribbons with spin-degenerate electronic structure and ne¬ 
glected spin-flip scattering, G splits equally, G|| = G/2 and 
G|| = G/2 [15]. In the case of magnetic ribbons, and i- 
current do not match anymore; unpaired spins in the ribbon, 
e.g. at defects or zigzag edges, can cause characteristic differ¬ 
ences between G^^ and Gn [16]. 

The spin-flip conductance coefficients G|| and G^^ are non¬ 
vanishing in the presence of (i) SOI [17] or (ii) exchange in¬ 
teraction with local spins in the device [18]. To include SOI 
in our DFT formalism, we employ an all-electron SOI mod¬ 
ule [19], the exchange-interaction is dealt with on the level of 
spin DFT [20]; details see Sec. II. 

As one would expect, we find a very small spin-flip con¬ 
ductance Go-o- in clean armchair GNRs (AGNRs) due to the 
very weak SOI and the absence of local impurity spins, see 
Sec. Ill A. In contrast, the spin-flip conductance is massively 
enhanced in the presence of adsorbates, see Sec. IllB - HID. 
For instance, our results indicate that the spin-flip probability 
associated with a single hydrogen adatom can be compara¬ 
ble to the spin-conserving one. This high spin-flip probability 
is rationalized by employing a simplistic tight-binding (toy) 
model. Our first-principles results are qualitatively similar to 
analytical results by earlier authors. Ref. [7], who employ a 
model calculation that is valid in the highly dilute limit. 


II. METHOD 

In our calculations, we are employing an extension of the 
AITRANSS platform, our DFT-based transport simulation 
tool [21, 22]. The spin-dependent conductance is obtained 
as follows: We extract the Kohn-Sham (KS) Hamiltonian 
H = [(77||, 77||), (Hi'i, Hii)], a 2x2 block-matrix in spin space, 
from a DFT calculation [23] including all-electron SOI [19] 
for a finite-size hydrogen-terminated graphene nanoribbon 
with horizontal armchair edges, see Fig. 1. Subsequently, we 
obtain the (retarded) single particle KS-Green’s function G of 
a finite-size strip in the presence of the left and right contacts 
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by standard recursive Green’s function techniques [14]: 



The spin-diagonal self-energies = [(S^, 0), (0, Sp] with 

reflect the presence of the leads [24] . They are treated 
with a closed-shell electronic structure and a vanishing SOI so 
that spin is a good quantum number in the leads [25]. Go rep¬ 
resents the bare KS-Green’s function of the device region, see 
Fig. 1 (a). We compute the spin-dependent zero-bias conduc¬ 
tance Gcrcr'iE) at a given chemical potential E [see Fig. 1 (b)] 
in a Landauer-Biittiker approach [26] : 

G^AE) = j Tr[f^ f« G^,^], (2) 

withf“=/(Z“-(S“)t). 

For spin quantization n = (sin 6 cos (f, sin 6 sin (f, cos 6) de¬ 
viating from default z-direction, we rotate the Green’s func¬ 
tion in spin space by the unitary transform U : 


G« = VGA, U = 


t 


f 


COS 


sm: 




sin I 


cos i 


(3) 


Due to the closed-shell electronic structure of the leads, fj re¬ 
mains unchanged by a unitary transform and the conductance 
with respect to an arbitrary spin quantization axis n is given 

by 2 

(£) = \ Tr[f^ (G« )t tl. G« 1 ■ (4) 


The formalism outlined here is well established [17]. 


III. RESULTS 

In this section, we present simulation results of the spin- 
dependent conductance Gera' of a clean and hydrogenated AG- 
NRs calculated according to Eqs. (2) and (4). For computa¬ 
tional details, we refer to App. A. 


A. Clean ribbon 



left lead 


right lead 



E-E^ (eV) 

FIG. 1. (a) Structure of a clean AGNRll device. In the blue marked 
device region, SOI is present, (b) Between the two leads of (a), an in¬ 
finitesimal voltage dV is applied. As response, an electronic current 
dl flows, which is splitted into four components dlo-A • They account 
for an electronic current being injected with spin cr and measured af¬ 
ter passing the device region with spin cr'. The spin-dependent con¬ 
ductance (E) at a given chemical potential E is then defined as 
ratio of dE^ and dV. The reference energy Fp is the chemical poten¬ 
tial of the isolated, charge-neutral device, (c) Corresponding conduc¬ 
tance of the pristine ribbon (a) according to Eq. (2). The spin quan¬ 
tization axis is chosen in positive z-direction. The spin-conserving 
conductance does not vanish completely inside the bandgap (but is 
reduced by seven orders of magnitude) due to the limited horizontal 
size of the lead system, so that the one-dimensional band structure of 
the ribbon cannot be formed with ultra-high precision. 


We calculate the conductance coefficients Go-o-'iE) of an 
AGNR with Ac = 11 transverse carbon atoms (AGNRll), 
see Fig. 1 (a) for the molecular structure of the device and 
Fig. 1 (b) for a sketch of the electronic structure and the trans¬ 
port setup. The result for the spin-dependent conductance ac¬ 
cording to Eq. (2) is shown in Fig. 1 (c) with a spin quantiza¬ 
tion axis in -rz-direction. For the spin-conserving conductance 
coefficients G^-^-, we find a step function with Go-AE)lie^ Ih) 
simply counting the energy bands intersecting with a given 
energy E [27]. Gcrcr is hardly affected by the SOI. Most im¬ 
portantly, the spin-flip conductance is found to be very small 
with an upper bound of 10“^^ jh due to the very weak SOI 
as expected. [28] Due to vertical mirror symmetry, and 
bands of pristine AGNRs are degenerate [29] and therefore 
G^^{E) = GiiiE) and G^i{E) = G^^iE) for z-quantization. 


B. Ribbon with a single hydrogen adatom 

We continue with a ribbon containing a single hydrogen 
adatom, see Fig. 2 (a) for the molecular geometry. Hydrogen 
forms a chemical bond with the underlying carbon atom re¬ 
sulting in an sp^ hybridization. The four nearest carbon atoms 
were structurally relaxed in order to catch the massive en¬ 
hancement of SOI due to the lattice distortion [30]. 

First, we comment on the finite-size DFT calculation 
of the ribbon including SOI as sketched in Fig. 2(b): 
The computed magnetization of the finite-size ribbon is 
(S) = (-0.16, -0.01, \A9Yfi [31], so there are approximately 
three unpaired electrons in the ribbon distributed at the zigzag 
edges [12] and near the impurity, see Fig. 2 (b). [32] 
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FIG. 2. (a) Structure of the AGNRll device with one single hydro¬ 
gen adsorbate, (b) Finite-size ribbon for the underlying DFT calcula¬ 
tion including SOL All atoms in the blue box belong to the device re¬ 
gion. The red surface denotes the isosurface of the spin density calcu¬ 
lated by the DFT including SOI where we observe the famous zigzag 
edge magnetism [12]. The total spin is calculated as (S) = f S(r)d^r. 
(c)-(f) Conductance of device region (a) with/without SOI and 
quantization axis according to inset. The reference energy, Ep is the 
chemical potential of the isolated, charge-neutral lead. The average 
spin-flip conductance Ga-a is computed as arithmetic mean of G^^ 
and G^ in the energy interval [- 1.5 eV, 1.5 eV]. 


Our simulation results for the spin-dependent conduc¬ 
tance Gcrcr' according to Eqs. (2) and (4) are displayed in 
Fig. 2 (c) - (f). First, we focus on the situation where the spins 
of sample and incoming electrons are (very nearly) aligned. 
The corresponding spin-conserving conductance Gera- is dis¬ 
played in Fig. 2 (c) with SOI and (e) without SOI. The values 
for G'i'^(E) and G^(£’) in Fig. 2 (c) and (e) deviate by less than 
I0~^e^lh, only, and we observe that SOI hardly influences 
Go-o-- In particular, the broad antiresonance [33, 34] indicat¬ 
ing the quasilocalized state (zero mode) that accompanies the 


isolated adatom remains clearly visible also in the presence of 
SOI. In Fig. 2 (c), the spin-flip conductance is seen to be very 
small, Gcrd-(E) ~ 10“^ e^/h, while it vanishes in Fig. 2 (e) due 
to the absence of SOI. 

Next, we consider in Fig. 2 (d) incoming electrons with spin 
polarization along x-axis, i.e. perpendicular to the magnetiza¬ 
tion axis of local spins in the ribbon. Here, the spin-flip 
conductance increases strongly reaching values larger than 
0.05 e^/h. Near the bandgap, it is even exceeding the spin- 
conserving conductance. In order to emphasize that the large 
spin-flip conductance is due to the exchange interaction and 
not related to SOI, we repeat the same calculation without 
SOI. The result is shown in Fig. 2 (f) and indeed it is indistin¬ 
guishable from Fig. 2 (d). In a nutshell, the exchange-driven 
spin flip is understood as follows [35]: Say, the flxed impurity 
spin points into the z-direction, i.e. the exchange interaction 
turns into S ^ ~ o-^ with the Pauli matrix acting 

on incoming conduction electrons. As a consequence, the ef¬ 
fective single-particle Hamiltonian no longer commutes with 
^cond’ of incoming particles no longer is 

conserved, if it happens to exhibit a component perpendicu¬ 
lar to the impurity spin. Therefore, spin-flips become possible 
with a probability (G^^ + G^^)/ Xo-o-' Gcrcr' that can reach or¬ 
der unity for a non-collinear spin passing a single hydrogen 
adatom. In App. B, we explain how our results are rational¬ 
ized employing a simple toy model. Our overall flndings are 
consistent with Ref. [7] that has employed a model calcula¬ 
tion. 


C. Ribbon with two hydrogen adatoms 

We calculate the spin-dependent conductance of a ribbon 
with two neighboring hydrogen adatoms, see Fig. 3 (a). This 
double-hydrogen defect is non-magnetic and we would like 
to confront it with the case of an isolated hydrogen adatom. 
In Fig. 3 (a), the non-magnetic character of this impurity is 
evident: The spin density near the impurity is smaller than 
the isovalue and so it cannot be resolved anymore. There are 
two reasons for the absence of magnetism with this defect, 
(i) There is no imbalance between the graphene sublattices, 
i.e. one impurity on each sublattice. This also implies that the 
number of electrons remains even, (ii) Because of the vicinity 
of the border, orbital degeneracies are lifted, so that a closed- 
shell ground state is favored. 

The computed conductance is shown in Fig. 3 (b) for the 
collinear case and in Fig. 3 (c) for the non-collinear one. In¬ 
deed, the spin-flip conductance for non-collinear transport is 
two orders of magnitude smaller as compared to the previ¬ 
ous case with a single hydrogen adatom. Note, that again, 
the spin-flip conductance for electrons polarized perpendicu¬ 
lar to the sample magnetization exceeds the collinear one - 
probably due to weak residual magnetism with spin density 
|S(r)|< 0.005 a.u.. 

We conclude that the interactions between hydrogen 
adatoms can be relevant, if they come sufficiently close. In 
our example, the spin-flip conductance of the ribbon with two 
hydrogen adatoms is far less than twice the value of the iso- 
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FIG. 3. (a) Finite-size ribbon with two hydrogen adatoms for the 
underlying DFT calculation including SOL The atomic structure is 
relaxed including the six surrounding carbon atoms. All atoms in the 
blue box belong to the device region, (b) Spin-dependent conduc¬ 
tance for z-polarized electrons and (c) for v-polarized electrons. 


FIG. 4. (a) Finite-size ribbon with 12 hydrogen adatoms for the un¬ 
derlying DFT calculation including SOI. The total number of un¬ 
paired electrons is Aspm = 2//z |(S)| = 1.75. All atoms in the blue box 
belong to the device region, (b) Spin-dependent conductance for 
electrons polarized along and (c) for x-polarization. 


lated adatom. This indicates a breakdown of Matthiesen’s rule 
in the limit of higher concentrations. 

D. Massively hydrogenated ribbon 

In Fig. 4, we show the spin-dependent conductance for 
a ribbon with 12 adsorbed hydrogen atoms which are dis¬ 
tributed randomly on the ribbon. The part of the structure 
that carries 44 carbon atoms and 12 adsorbed hydrogen atoms 
was structurally relaxed. The transmission function displayed 
in Fig. 4 (b) and (c) is seen to carry strong mesoscopic fluc¬ 
tuations that reflect many quasilocalized states near the Fermi 
energy. In this situation, the spin-flip conductance can reach 
the same order of magnitude as the spin-conserving one in 
sizable energy window. 

Notice, that even for incoming electrons with collinear spin 
along the direction of (S), the spin-flip conductance is very 
large. We interpret this effect as an indication that the direc¬ 
tion of the local spin density is fluctuating in space. 

E. Local current density in massively hydrogenated ribbon 

In Fig. 5, we show the local current density response for 
the massively hydrogenated AGNRll of Sec. HID based on 
an open-shell DFT calculation including SOI (see App. C for 
method details; see Fig. 4 (a) for the atomic structure). 

The current exhibits strong mesoscopic fluctuations cover¬ 
ing three orders of magnitude. They are related to vortices 
which exceed the average current by over one order of magni¬ 
tude. Such current vortices go along with orbital magnetism 
which are also potentially relevant for spin relaxation [8]. 



FIG. 5. Local current density response (integrated over the out- 
of-plane direetion; normalized to the eonduetance per width) in the 
hydrogenated AGNRl 1 shown in Fig. 4 (a) (blue box only). The cur¬ 
rent density exhibits strong mesoscopic fluctuations that reflect in a 
logarithmic color scale covering three decades. Some interesting cur¬ 
rent paths are drawn into the picture for illustration: Local current 
vortices exceeding the spatial average current by one order of mag¬ 
nitude (see dark red regions). Plot shows current amplitude (color), 
current direetion (arrows), carbon atoms (grey crosses) and hydrogen 
atoms (red crosses). 


We compare this finding with our recent work [8], in which 
we simulated the local current density for a larger ribbon 
[Nc = ^l transverse carbon atoms (AGNR41)] but enforced 
a closed-shell electronic structure without SOI to reduce the 
computational effort [36]. The characteristic features, i.e., 
current vortices and broad fluctuations, are the same in both 
cases; they emerge from the full spin treatment including SOI 
as well as from the spin restricted treatment. Therefore, we 
confirm that their appearance is a fundamental property of the 
scattering states in defected graphene flakes. 
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IV. CONCLUSION 

In conclusion, we calculate the spin-dependent zero-bias 
conductance Gg-o-' in armchair graphene nanoribbons (AG- 
NRs) with hydrogen adsorbates employing a DFT-based ab 
initio transport formalism including spin-orbit interaction 
(SOI). We find that a narrow AGNR decorated with a sin¬ 
gle hydrogen adatom exhibits a spin-flip conductance Gg-o- 
that is highly anisotropic. In the case of collinear conduct¬ 
ing and local impurity spins, the spin-flip conductance is due 
to SOI and it is very small, at most In contrast, 

in the non-collinear situation, we observe a spin-flip conduc¬ 
tance between 10“^ e^/h and 10“^ e^/h, that can even exceed 
the spin-conserving conductance in some energy range. We 
explain this effect by an exchange-mediated spin-flip mecha¬ 
nism masking any potential spin-flip effect by spin-orbit in¬ 
teraction. Our calculations suggest that the spin-dependent 
conductance becomes isotropic again, if the concentration of 
adatoms is not too small. In this case the exchange mediated 
spin-flip scattering is always strong. 
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APPENDIX: TOY MODEL, COMPUTATION OF 
CONDUCTANCE AND CURRENT DENSITY 

In App. A, we present details on the computation of the 
spin-dependent conductance. In App. B, we introduce a toy 
model to rationalize the shape and the order of magnitude of 
the spin-dependent conductance of the ribbon with a single 
hydrogen adatom. Appendix C provides information on the 
formalism to calculate the current density shown in Fig. 5. 


Appendix A: Computation of the spin-dependent conductance 

1. Method details: Expansion of operators in real-space basis 
functions and partitioning in device and contact region 

In this section, we describe the expansion of all operators 
appearing in Eqs. (1) - (4) in basis functions of the underlying 
DFT calculation and how we partition the appearing matri¬ 


ces in blocks belonging to the device region and the contact 
regions to the reservoirs. 

From a non-periodic open-shell DFT calculation includ¬ 
ing all-electron spin-orbit coupling (SOI) [19], we obtain the 
Kohn-Sham (KS) matrix H of a finite-size ribbon [e.g. see 
Fig. 2 (b)] with matrix elements 

= J dh. (Al) 

denote the basis functions constructed from the un¬ 
derlying DFT calculation [37] which are real-valued, atom- 
centered and orthogonalized (via Lowdin orthogonaliza- 
tion [38]). H is of size 2Nx2N with non-zero entries in the 
off-diagonal spin-blocks due to SOI. 

Subsequently, we cut off every element from H, if 
and/or ^ is centered on an atom outside the device re¬ 
gion, the latter indicated by the blue boxes in the previous fig¬ 
ures. Employing the resulting truncated 2Ax2A-device-KS 
matrix Hdevice {N <N), we calculate the 2Ax2A(-matrix rep¬ 
resentation of the device-Green’s function in presence of the 
left and the right reservoirs as 

G{E) = [^1 

2N Hdevice . (A2) 

The 2Ax2A(-self-energy matrices Il^(£’) are computed by a 
separate DFT calculation of a clean, finite-size and closed- 
shell treated AGNRll ribbon, see Ref. [22] for details. The 
self-energy matrices TE{E) are spin block-diagonal because 
of the closed-shell treatment of the reservoirs. Their non¬ 
vanishing entries correspond to the contact regions (marked 
by magenta boxes in Fig. 6) which - throughout this paper - 
were chosen as the outermost left and the outermost right col¬ 
umn of carbon rings inside the (blue marked) device-region 
(i. e. 22 carbon atoms belong to each of the left and right con¬ 
tact region). 

Then, we partition G(£') into four AxA-block matrices 
Gcrcr'{E) and Il^(£’) into two AxA-block matrices 
[shorthand notation: Il^(£’) with Il^(£’) = due to the 

closed-shell electronic structure of the reservoirs]. Finally, 
the spin-dependent conductance Gcrcr'iE) is computed as an 
orbital trace of a product of AxA-matrices, 

G^^.(^)= -Tr[r^(£)(G^^.(^))^r^,(^)G^v(^)] , (A3) 
where r“(£) = /(!:“(£) - (2:“(£))^). 

2. Validiation of the partitioning in device and contact region 

In this section, we show that the spin-dependent conduc¬ 
tance is only weakly dependent on the precise (numerical) 
partitioning of the finite-size ribbon into device and contact 
region. As an example, we consider a ribbon with a single 
hydrogen adatom [as in Fig. 2 (a)]. As finite-size input geom¬ 
etry for the SOI-DFT calculation, we choose a long finite-size 
ribbon, see Fig. 6 (a)/(b), but we define two different device 
regions, as sketched by the blue boxes. The contact regions 


6 



o device region O contact region 
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As discussed in the body of this paper, the high spin-flip 
conductance is due to exchange with the magnetic moment of 
the device region. As can be seen in Fig. 2 (b), the unpaired 
electrons are localized near the zigzag edges and the hydrogen 
impurity. For a calculation of Gcrcr'(E) without any finite-size 
artifacts, the full local moment near the hydrogen adsorbate 
has to be included in the device region, but no local moments 
originating from the zigzag edges. We conclude by compar¬ 
ing Gcrd-(E) in Fig. 6 (c) and (d), that local moments of the 
zigzag edge states are included in the large device region of 
Fig. 6 (a) and/or that not the whole local moment caused by 
the hydrogen adsorbate is contained in the small device re¬ 
gion of Fig. 6 (b). However, the difference between Gcrd-(E) 
of Fig. 6 (c) and (d) is not serious as mentioned above. Con¬ 
sequences of changing the contact region are investigated in 
Ref. [22]. 


Appendix B: Toy model for spin flips due to local exchange 

To explain the shape and order of magnitude of Gg-o-' (E) of 
Fig. 2, consider a toy model that consists of the Hamiltonian 
of the lead and a single site 

Ho = - to ^ ^ (^I+l,cr^-^’Cr + c\crCx+\,o^ (^1) 

^ Sg-dl-dcr + (B 2 ) 

and a coupling term 


FIG. 6. (a) Input geometry for the SOI-DFT calculation of an 

AGNRl 1 with a single adsorbed hydrogen atom. The computed total 
magnetic moment of the ribbon is (S) = (- 0.71, - 0.48, - QSlYfi, so 
that there are 2.7 unpaired electrons in the ribbon. For the transport 
calculation, we partition the ribbon into a device region (blue box) 
and two contact regions (magenta boxes), (b) Identical geometry and 
identical DFT calculation as in (a), but different partitioning in de¬ 
vice and contact regions. (c)/(d) Result for the spin-dependent con¬ 
ductance according to Eq. (A3) for partitioning (a) [shown in (c)] and 
for partitioning (b) [shown in (d)] for z-quantization so that e(S) ^ 
enabling the efficient exchange-mediated spin-flip mechanism. 


with non-vanishing self-energy are chosen as outermost row 
of carbon rings in the device region as indicated by the ma¬ 
genta boxes. 

The computed conductance Gg-o-' {E) for the large device re¬ 
gion of Fig. 6 (a) is shown in Fig. 6 (c) and for the small device 
region of Fig. 6 (b) in Fig. 6 (d). The spin-conserving conduc¬ 
tance Gcrg-{E) in (c) and (d) are indistinguishable from each 
other, while the spin-flip conductance Gcrd-{E) in (c) and (d) 
exhibit the same shape. They are only deviating quantitatively 
from each other in energy intervals, where Gcrd-{E) is anyway 
small. Additionally, the order of magnitude of Gcrg-{E) agrees 
to the one in Fig. 2 (d)/(f). 


t> = fLC X + 

with Co- = co,o-- The whole model is sketched in Fig. 7 (a). 
It resembles a hydrogen adsorbate that splits off a resonant 
level [33, 39] near the charge-neutral point from the conduc¬ 
tion band continuum for both spin channels. The exchange 
interaction with the local spin lifts the degeneracy between 
the localized states with differing spin: In the case 

of a non-vanishing overlap matrix element the associated 
quasilocalized state contributes a separate conductance chan¬ 
nel that interferes destructively with the residual ones [40]. 
The SOI is modelled to be only present at the resonant site 
by the parameter A. The model Eqs. (B1)-(B3) account for 
collinear conducting and local spins reflecting in the vanish¬ 
ing overlap between cr band and a quasilocalized state. The 
spin-dependent conductance of the collinear-spin model of 
Fig. 7 (a) is depicted in (c). We observe the well-known con¬ 
ductance dips at the resonant energies and si [33] and a 
small spin-flip conductance in the order of (T/O^ due to the 
small SOI parameter A. 

When rotating the spin quantization of incoming elec¬ 
trons by 7r/2 from yi to -r/-, |+) = (IT) + U))/V2 and 
l-) = (IT>-li))/V2, an overlap of both quasilocalized states 
to both bands is formed. To account for this we generalize our 
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(a) it)^ |t>^ (b) 1+)^ 1+) 



FIG. 7. (a)/(b) Tight-binding model to rationalize the conductance 
curves of the AGNR with a single hydrogen adsorbate including two 
conducting bands, a localized state with spin t and J, and a SOI 
between both localized states. The matrix elements between the 
quantum wire and the localized states are (a) ^lc in the case of in¬ 
coming electrons being spin-polarized along t H and (b) ^lc/V 2 for 
+! —polarization. Here, one hopping term is negative, see Eq. (B6). 
(c)/(d) Spin-dependent conductance for model parameters for the 
AGNR with hydrogen adatom and (c) t/i- and (d) +!—quantization 
of incoming electrons. The SOI strength X is estimated by the SOI 
strength of 2.5 meV for a hydrogen adatom on graphene [30] in rela¬ 
tion to the graphene hopping t = 2.8 eV. 


model 

Ho = - to ^ ^ (B4) 

fi=+,- X 

Ho= Scrdl-dcr + /T (B5) 

and the coupling 

y = -^ (ij (c+ + c_) + dl (c+ - c_) + h.c.) (B6) 

where c+ = + c^)/V2 and c_ = - c^)/V2 [model in 

Fig. 7 (b)]. It now imitates the situation where the local mo¬ 


ment and the conduction band spin are not collinear. Then, 
even in the absence of SOI, A = 0, the conduction band spin 
is not conserved and we expect a large spin-flip rate. Indeed, 
seen in Fig. 7 (d), the spin-flip conductance increases by four 
orders of magnitude compared to (c) even exceeding the spin- 
conserving conductance near the antiresonance. The results of 
the model agree well with the conductance of the ribbon with 
a single hydrogen adatom in Fig. 2: In particular, it also repro¬ 
duces the quasilocalized state also seen in or near the bandgap 
of the ribbon. 

We mention that in the non-magnetic case, £'^=£i, the 
efficient spin-flip process for -r/—polarization is suppressed. 
The reason is that in this case the single site effectively acts 
like a two-path interferometer that support perfect destructive 
interference for the two tunneling paths between the -f and 
the - band. 

Appendix C: Computation of the local current density 

For completeness, we summarize the calculation of the 
local current density response j(r, E) starting from the full 
spin-dependent retarded Green’s function G(E) [see Eq. (1) 
for its spin components]. The retarded Green’s function al¬ 
lows to calculate the non-equilibrium Keldysh Green’s func¬ 
tion G^(E): 

G< = /G[/LfL + , (Cl) 

with f“ = i(t" - (2“)^). The occupation numbers /“ of the 
leads reduce to step functions at zero (or low) temperature. 
Inside the voltage window, we assume an occupied left lead 
and an empty right lead, i.e., = 1, = 0, so that reduces 

to 

G^(E) = iG(E)tHE)G\E). (C2) 

The Keldysh Green’s function is transformed to real-space 
representation using the basis functions of the underlying 
DFT calculation: G^^Xr,r', E) = {ro-\G^(E)\r'o-'). The cur¬ 
rent density (per energy) is then expressed as 

j(r, £) = 1A y iim(Vr, - Vr)Gy(r, r',E). (C3) 

271 2m r'^r 

cr 

The factor 27i reflects an inverse Fourier transform. 
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